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Abstract Observations of H2 spectra in the line-of-sight of distant quasars may 
reveal a variation of the proton-electron mass ratio fi = m p /m e at high redshift, 
typically for z > 2. Currently four high-quality systems (Q0347 - 383, Q0405 - 443, 
Q0528 — 250 and J2123 — 005) have been analyzed returning a constraint Aji/jj, < 
1 x 10~ 5 . We present data and a /z-variation analysis of another system, Q2348 — 
Oil at redshift Zabs = 2.42, delivering A\ij\i = (-1.5 ± 1.6) x 10~ 5 . In addition 
to observational data the status of the laboratory measurements is reviewed. The 
future possibilities of deriving a competitive constraint on Ajj./ji from the known 
high-redshift H2 absorbers is investigated, resulting in the identification of a number 
of potentially useful systems for detecting jU-variation. 
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1 Cosmological ji variation and the laboratory database 

The proton-electron mass ratio ji = m p /m e is one of the constants suitable for tar- 
geting a possible variation of a physical constant on a cosmological time scale. Al- 
though ji is a dimensionless constant, lending itself for an observational approach, 
it is not a. fundamental constant in a strict sense, like the fine-structure constant a, 
because jj. involves all the binding physics of the quarks inside the proton. How- 
ever, its close connection to ;« c /AgcD makes it the parameter determining the mea- 
sure of the strong force with respect to the electroweak sector, and is therewith of 
fundamental importance. Variation of fi can be straightforwardly derived from a 
comparison between spectral lines observed at vastly differing redshifts. Values of 
wavelengths at high redshift XI are compared to laboratory wavelengths X' Q for as 
many as possible H2 lines [ 1| using sensitivity coefficients for each individual spec- 
tral line Ki = dlnXi/dlnji. These Kj are calculable to the 1% accuracy level 0. 

The laboratory wavelengths of the spectral lines in the Lyman and Werner bands 
of H2 have been investigated at ever improved accuracy using tunable vacuum ul- 
traviolet (VUV) lasers since over a decade [3, 4j [3J [6J [7] E ], also including spectra 
of HD (9] [TO]. In addition the novel Fourier-transform spectrometer in the VUV 
range at the Soleil synchrotron was used to perform direct absorption spectra focus- 
ing on HD to cover the entire range of Lyman and Werner bands for this deuterated 
molecule ifTTIl . The most accurate wavelength positions are obtained from a third in- 
direct spectroscopic method based on the highly accurate determination of level en- 
ergies for certain anchor levels fl2l . where the determination of Lyman and Werner 
wavelengths results from an independent measurement of spacings between excited 
states [13]. The full set of recommended wavelengths is reported in Bailly et al. Ifl4ll . 
Another compilation of these data, both for H2 and HD is given in the appendix of 
Malec et al. 03), also including a listing of Kj sensitivity coefficients (2l[TT], oscil- 
lator strengths of the absorption lines lfT6l . and of damping coefficients ifTTl . 

A quasar absorption spectrum can be treated in different ways to detect a possible 
Ajj./ji. A straightforward method is to derive line positions by fitting H2 resonances, 
and thereupon compare the results with the laboratory wavelengths; by deriving re- 
duced redshifts = (z,- — za)/ (1 + Za) with z; the redshift determined for each line 
(Zi = (A,' /A.?) — 1) and za the redshift of the absorber, such a method gives a graphi- 
cal insight into a possible variation of jj, J3] |7] . In contrast to this line-by-line fitting 
method a comprehensive fitting method has advantages. There, an H2 spectrum is 
generated, based on the available molecular parameters (wavelengths lfT4l[T5l . os- 
cillator strengths lfT6l . and damping parameters [ 17 1). The value for Aji/fi is then 
a single fitting parameter comparing the observed spectrum with a model spectrum; 
this is accomplished via least squares fitting with the VPFIT software ifTSl . Each 
line is convolved with the instrumental profile (pertaining to the spectral resolving 
power R — X/AX of the spectrometer), a natural linewidth F, and a Doppler pa- 
rameter b. The main advantage of the method is that a composite velocity structure 
can be addressed. The column densities N(J) for each rotational state can be tied 
between velocity components or left independent. Further details and advantages of 
the comprehensive fitting method are discussed in Refs. lfT5l[T9ll . By this means the 
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Fig. 1 Current result on cosmological /l-variation based on the five available high-quality H? ab- 
sorption systems as a function of redshift: blue points relate to Q0347 — 383, Q0405 — 443 and 
Q0528 - 250 1 19], green is the result from J2123 - 005 Q3), while the black points relate to the 
ammonia data 1201 1211 for redshifts z < 1 . The red hexagonal point represents the result of the 
present analysis of Q2348 — 01 1. 



spectral information on overlapping lines and of saturated lines can be included. In 
the example of the analysis of the Keck-HIRES spectrum of the quasar J2123 — 005 
it is shown how the velocity structure can be investigated in much detail in terms of 
computed composite residuals which demonstrate a complex velocity structure that 
should be accounted for [ 15 1. These optimized comprehensive fitting methods were 
also implemented in the analyses of the three systems Q0347 — 383, Q0405 — 443 
and Q0528 - 250 by King et al. 09). 

The resulting constraints on A for the four systems mentioned are shown in 
Fig. [T] The figure also inculdes the result of the Q2348 — 01 1 analysis, details of 
which will be presented in the next section. The tighter constraints stemming from 
radio-astronomical observations of NH3 at redshifts z < 1 are shown as well 02011211 . 
The data show the present status of a possible variation of the proton-electron mass 
ratio to be A\ij[i < 1 x 10~ 5 at redshifts z = 2 - 3.5. 
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2 VLT-UVES observation of Q2348-011 



H2 absorption features associated with the Damped Lyman-a (DLA) system at z = 
2.42 along the sightline of Q2348 - 01 1 (SDSS J235057.88 -005209.8, z em - 3.02) 
were first reported by Petitjean et al. Il22l . Since then this system has been the subject 
of several investigations but none of them focused on jU-variation (e.g. Ledoux et 
al. Il23l . Noterdaeme et al. 11241 ). but rather on the molecular column densities, the 
metallicity and the physical conditions of the gas contained in this DLA-system. 
For an investigation of /^-variation new observations of well-calibrated spectra of 
Q2348 — 01 1 were recorded, and will be presented here. 




-150 -75 75 
Relative velocity (km s _1 ) 

Fig. 2 Lower panel: H2 velocity structure in the Q2348 — 011 DLA system at z = 2.42; upper 
panel: Fe II velocity structure of the same absorber. The positions of 7 different velocity features 
of H2 are indicated by the blue dashed lines, while the positions of 14 velocity components of 
Fe II are indicated by the red dashed lines. Note that the strongest velocity component in Fe II 
absorption profile (at km s~') does not have an associated FL component. The cleanest Fe II 
transition (1608.45 A) is blended with Si IV (1402.77 A) of an additional absorber at z = 2.93. The 
green ticks in the Fe II panel indicate positions of the overlapping Si IV lines. 



2.1 Data 



The absorption spectrum of the quasar Q2348 —011 was obtained with the Ul- 
traviolet and Visual Echelle Spectrograph (UVES) at the Very Large Telescope 
(VLT) of the European Southern Observatory (ESO) in Paranal, Chile. Observations 
were carried out in visitor mode on four consecutive nights (2007 August 18-21). 
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Q2348 — Oil is not a particularly bright object (R mag = 18.31) so long exposures 
were needed. An overall of 15 spectra were recorded with an exposure time between 
3600 and 4800 seconds each; they make up for a combined total of 19.25 hours of 
observation. The slit widths were set to 0."8 — 1 ."0 and 0."7 — 0."8 in, respectively, 
blue and red ranges. They were kept constant for all exposures irrespective of the 
noticeably varying seeing conditions, which were on average l."2. These settings 
yield a resolution of R ~ 51 700 (in the blue-UV range where all H2 lines are ob- 
served) and 55 300 (in the red). This corresponds to a resolution element of 2.46 and 
2.31 km s _1 in terms of velocity units, respectively. For optimum wavelength cali- 
bration the 15 science exposures were each attached with a ThAr calibration frame, 
obtained immediately after each science exposure. Special care was taken to use 
pre-selected Th-Ar reference lines following optimization procedures l25l . The fi- 
nal spectrum of Q2348 — 01 1 covers (vacuum-heliocentric) wavelengths 3572-9467 
A, with gaps at 4520-4621 and 7505-7665 A. Parts of the spectrum are displayed in 
Fig.0 



2.2 Analysis 

The wavelength range covered by the VLT-UVES spectrum provides 58 H2 transi- 
tions for rotational levels 7 = — 5. The H2 features fall between 3580 and 3860 
A for this absorber at z = 2.42. This part of the combined spectrum has a signal- 
to-noise ratio (S/N) of 25. Note that no associated HD spectral features have been 
detected at a significant level. All the transitions observed arise from the Lyman 
band. Molecular hydrogen is present in 7 velocity features (the two left-most are 
blended together but they count at least as two), so in principle this spectrum would 
provide a sample of 58 x 7 absorption lines. However, besides the DLA at z = 2.42 
(neutral content (logA^(H I) = 20.50 ±0.10), ||23l ), which contains molecular hy- 
drogen, there is an additional DLA at z = 2.62, which has a larger column density 
of neutral hydrogen (logA^(H I) = 21.30±0.08) l24]| ). but no signs of molecular 
hydrogen. The strong Lyman-j3 feature of the second DLA (z = 2.62) at 3710 A 
obscures 4 H2 transitions (24 velocity components) and damps some others (see 
Fig. [3]), while also the Lyman-y blends 4 H2 lines (28 velocity components). This 
DLA gives rise to a cutoff in the absorption spectrum at 3580 A, which makes this 
DLA an unfortunate coincidence along the sightline of Q2348 — 011. In addition 
there exist strong H I absorbers, almost at the column densities of sub-DLA's, at 
z — 2.73 and z = 2.93 obscuring several H2 lines. 

Another distinctive feature of this spectrum is self-blending of H2: since the ve- 
locity structure of H2 consists of numerous components spread over ~ 300 km s , 
nearby lying lines exhibit mutual overlap of velocity components. For example, ve- 
locity components 4-7 of H2 L0R1 lie in the same region as velocity components 
1-3 of H2 L1P5. This kind of overlapping can be disentangled using the compre- 
hensive fitting method, if a sufficient number of non-overlapping lines is available. 



6 W. Ubachs, J. Bagdonaite, M.T. Murphy, R. Buning, and L. Kaper 




3720 3725 3730 3735 3740 3745 3750 



1 








3750 3755 3760 3765 3770 3775 3780 



Fig. 3 Hi absorption model superimposed on a portion of the Q2348 — 01 1 spectrum. The strong 
feature centered at 3710 A is Lyman-/3 of the DLA at z = 2.62. 

However, the non-overlapping lines are blended in addition with the absorbers of 
the Lyman-a forest, and the overall amount of clear features is strongly reduced. 

Some of the molecular hydrogen features are contaminated by metal lines asso- 
ciated with the two DLAs, namely, L2R0 with Si II (1020.69 A) at z = 2.62, and 
L3R1, L3P1, L1P2, L1R3, L0P2, L0R3 are blended with various Fe II lines asso- 
ciated with the z = 2.42 DLA. Further, L1R5 and L1P4 may be blended with an 
additional C IV system at z = 1 .44 but this is difficult to reconstruct unambiguously. 
Generally, the sightline towards the quasar Q2348 — 011 is rich in interesting ab- 
sorbers: besides the two DLAs there are other strong metal absorbers at z = 2.97, 
2.93, 2.72, 2.58 (C IV, Si IV, N V), and at z = 0.86, 0.77 (Mg II). The latter is 
remarkable in having numerous components spread over 750 km s . Noteworthy 
coincidence is that the cleanest Fe II transition (1608.45 A) associated with z = 2.42 
is blended with Si IV (1402.77 A) of the absorber at z = 2.93 (see Fig.^. Also note 
that the most prominent velocity component in metal transitions does not have an 
associated H2 component. 

For the final fit deriving Afi/jj, 35 regions with 36 different molecular hydrogen 
transitions have been selected. Each line is modeled as a Voigt profile. In the com- 
prehensive fitting method the parameters (column density N, Doppler parameter b, 
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and position z) are connected for different lines. Each /-level for H2 has a differ- 
ent ground state population, so all transitions from the same /-level have the same 
value for N, b, and z. It is also assumed that a given velocity component has the 
same z and b value in all /-levels; physically this can be interpreted that the molec- 
ular hydrogen of all states resides in the same absorbing cloud. The value for Ajj./jj. 
is a single fitting parameter for all molecular hydrogen transitions. The absorption 
model is optimized to fit data by using the VPFIT software. 

As analyzed by Malec et al. in the example of the J2123 — 005 quasar system, 
more than one velocity component can underlie a molecular absorption feature as 
distinguishable by inspection. Their presence can be confirmed or rejected rely- 
ing on goodness-of-fit and by inspecting composite residuals 1151 . As the two left- 
most features of H2 absorption are blended, more components could in principle be 
present close to or between the two blended features. A fit with 7+1 velocity com- 
ponents (VC) instead of 7, where the additional VC was included between the two 
left-most VCs, has improved the Xv P er degree-of-freedom by several hundredths, 
but the additional velocity component was rejected as unnecessary in transitions 
from / = 4 and / = 5 rotational levels. The 7+1 VC model is regarded as the fidu- 
cial result, because it gives the best match with the observed spectrum, i.e. it is 
statistically preferred as it returns the smallest Xv- A more detailed analysis will be 
presented in an upcoming paper (27). 



3 Results 

The fiducial absorption model led to the following result: (Apt/ pL)i+ivc — (— 1.5± 
1.6) x 10~ 5 . Table [I] provides molecular cloud properties with their statistical un- 
certainties as derived from the final fitting attempt when Aji/jj. was set as a free 
parameter. 

The analysis of H2 absorption in the Q2348— 01 1 spectrum yields a competitive 
result on Aji/jj,, but it is three times less tight compared to the constraints from pre- 
vious analyses (see Fig. [TJ. One of the major causes that led to the relatively large 
uncertainty is the low S/N of the spectrum. The second cause is that the spectrum is 
obscured by the neutral hydrogen features of the additional strong DLA at z = 2.62. 
The H2 transitions falling at the shortest wavelengths covered by the spectrum are 
especially useful in /x-variation analysis, since they exhibit the larger /f-coefficients 
(transitions in the Lyman band with higher vibration excitation). In the spectrum of 
Q2348 — 01 1 they are not detected due to the H I absorption produced by the addi- 
tional DLA. Although the total number of H2 transitions available in this spectrum 
is reduced, it is at least partially compensated by the presence of numerous velocity 
components associated with each transition. Altogether, each transition contributes 
a higher information content when it is imprinted in several velocity components. 
The complex absorption structure was successfully modelled by means of the com- 
prehensive fitting method. 
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4 Potential high-redshift H2 absorbers for future detection of 
ji -variation 

Detection of cosmological /x-variation, or putting strong constraints on the variation, 
should preferably be based on large numbers of high redshift H2 absorbing systems. 
For comparison the recent study on a possible spatial dipole for the fine structure 
constant a relies on the analysis of almost 300 absorbing systems [26|. However, for 
H 2 , up to now only five systems (Q0347 - 383, Q0405 -443, Q0528 - 250, J2123 - 
050 and Q2348 — 011) of sufficient quality are found, and the results obtained in 
comprehensive fitting procedures are summarized in Fig. [T] We note that alternative 
line-by-line fitting procedures on the Q0347 — 383 and Q0405 — 443 systems have 
recently been reported [28 29 30]. Taking the weighted mean of the H2 results of 
Fig.[T]yields Aji/ix = (2.9 ±2.6) x 10~ 6 for redshift ranges between z = 2.0 -3.1, 
which may be considered as the status by September 2010. The data provide firm 
evidence that /z-variation is below the 10~ 5 level for look-back times in the range 
of 10-12 Gyrs. 




20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 

RAJ2000 / degrees 



Fig. 4 Known quasars with magnitude V < 18 as obtained from |32|. 

A recent analysis of the J1337 + 315 spectrum by Srianand et al. [j 3 1 1 yielded a 
constraint only at the 10~ 4 level. This is due to the combination of low brightness of 
the quasar (R mag = 18.08) and low H2 content at the absorbing cloud (log 10 .iV(H2)= 
14. 1). The present analysis of the Q2348 — 01 1 H2 absorption system demonstrates, 
that even when the S/N of a spectrum is low, a competitive constraint on A ji/jj. can 
still be delivered if other circumstances are favorable. In case of Q2348 — 01 1 the 
low S/N is compensated by the presence of at least 7 sufficiently strong velocity 
components. However, in this regard the Q2348 — 01 1 absorbing system is rather an 
exception than a rule. In other known DLAs H2 absorption is usually detected in no 
more than three visually distinguishable spectral features. 

This raises the question how many of the thousands of quasar systems on the 
sky (see Fig.|4]and |32|), of which there are some 1000 identified as DLAs could 
produce a competitive result, i.e. at the level 8(Ajx/jj.)<5x 10~ 6 . Based on our ex- 
perience in analyzing spectra, we postulate that high quality H2 absorption systems 
should obey the following conditions: 
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• An flux of V mag < 18 of the quasar background source for providing a S/N of 
about 50 after a reasonable amount of data collection (15-20 h). 

• A column density of iV(H2) preferably in the range between 10 14 and 10 18 cm 
yielding sufficient but not too saturated H2 absorption. 

• An H2 absorption system at redshift z > 2; only then a sufficient number of lines 
will shift into the atmospherically transparent window and the throughput win- 
dow of Keck-HIRES and VLT-UVES; preferably some Werner lines should be 
observed for dealing with systematics [ 15 1. 

• No occurrence of additional DLAs and only a limited amount of metal absorption 
systems along the same line of sight; the Q2348 — 01 1 absorber is an example of 
an unfavorable situation with an additional strong DLA and a few very strong H 
I absorbers, as well as several metal bearing systems. 

• In principle, a firm constraint on the velocity structure of the H2 absorber is only 
possible with the comprehensive fitting method. The present example of Q2348 — 
Oil demonstrates the applicability of this method in case of a complex structure 
with 7 velocity components. For line-by-line fitting methods the H2 spectrum 
should preferably exhibit a single absorber; that a, e.g. two-component model, is 
statistically worse than a single component model should be demonstrated. 

We have inspected the existing literature on high-redshift H2 absorbers in order to 
identify possibly useful systems that may be targeted in future for /i-variation. In 
Table [2] all systems in which H2 has been detected so far are listed. Relevant details 
for all systems, such as absorption and emission redshifts, their position on the sky, 
the obtained column densities for H2, HD and H I, and the magnitude are listed as 
well. Inspection of these systems leads to a quality assessment of the known H2 
absorber systems. Some of the systems can be discarded based on their low redshift 
(Q0551 - 366, Q2318 - 1 1 1, Q0013 - 004 and Q1331 + 170), on their low H 2 col- 
umn density (Q2343 + 125 and Q1337 + 315), too low brightness, or combinations 
thereof. It should be noted that this assessment is linked to the status of present day 
technology. Future larger telescopes, or space-based telescopes equipped with high 
resolution spectrometers {R > 50000) may alter this view. 

Based on the considerations given we have identified five quasar systems that 
have the potential to achieve a constraint at the 5(4/x//i)<5x 10~ 6 level. These 
systems have been investigated to some extent, but larger numbers of exposures 
yielding better S/N, and ThAr attached calibration frames will be needed to reach 
the desired accuracy. Of the five systems, four have been observed with VLT-UVES, 
while one (FJ0812 + 320) has been studied with Keck-HIRES. 

• Q0642 506 with redshift z a bs = 2.66 provides a good number of H2 lines with 
a high column density giving rise to saturation of only the strongest lines. The 
occurrence of only a single velocity components renders the spectrum relatively 
simple, also accessible for a line-by-line treatment. Spectra are shown in l33l . 
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• HE0027 184 is of similar good quality. It has a favorable magnitude and a 
molecular column density which allows for observation of a large number of 
lines to be observed. Indeed in the archived spectra with data observed in 2004 
many lines are found, possibly even up to J = 6 [34|. 

• FJ0812+320 has a strong saturated H2 component, which makes it less ideal. 
The system has been studied in particular for its abundant presence of the widest 
variety of metal lines [35 | and of HD l36l [37l . The system is nevertheless of 
importance because it is the only system of reasonable quality at northern decli- 
nation, except for the special z — 4.22 system discussed below. 

• J1237+064 at z a bs = 2.69 shares some of the disadvantages of the previous sys- 
tem, in particular the very high column density giving rise to saturation. However, 
the velocity structure exhibits two weaker components, such that high-/ lines can 
be obtained in the strong components, and the low-/ lines in the weak compo- 



Table 2 Listing of known high redshift H2 absorption systems with some relevant parameters. 
Bessel R photographic magnitude taken from the SuperCOSMOS Sky Survey (40). The five high- 
quality systems analyzed so far are offset. The column densities 7V(H2), Af(HD) and N(ti I) are 
given on a log 10 scale in [cm -2 ]. 
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+32:08:08.6 


19.9 


15.4 


21.4 


17.88 


HD ED 


Q0841+129 


2.37 


2.48 


08:44:24.24 


+ 12:45:46.5 


14.5 




20.6 


17.64 


ED 


Q1232+082 


2.34 


2.57 


12:34:37.58 


+07:58:43.6 


19.7 


15.5 


20.9 


18.40 


ED EE) 


J1237+064 


2.69 


2.78 


12:37:14.60 


+06:47:59.5 


I9.2 d 


14.5 


20.0 


18.21 


(38) 


Q1331+17<K 


1.78 


1.78 


13:33:35.81 


+ 16:49:03.7 


19.7 


14.8 


21.2 


16.26 


EHUD 


Q1337+315 


3.17 


3.17 


13:37:24.69 


+31:52:54.6 


14.1 




21.4 


18.08 


ED 


Q1439+113 


2.42 


2.58 


14:39:12.04 


+ 11:17:40.5 


19.4 


14.9 


20.1 


18.07 


ED 


Q1443+272 


4.22 


4.42 


14:43:31.18 


+27:24:36.4 


18.3 




21.0 


18.81 


|39| 


Q1444+014 


2.08 


2.21 


14:46:53.04 


+01:13:56.0 


18.3 




20.1 


18.10 


|48| 


Q2318-111 


1.99 


2.56 


23:21:28.69 


-10:51:22.5 


15.5 




20.7 


17.67 


ED 


Q2343+125 


2.43 


2.52 


23:46:25.42 


+ 12:47:43.9 


13.7 




20.4 


20.18 


E3|49) 



" A recent re-analysis has shown that this system contains HD 1501 . 

* Multiple DLAs in this system. 

Studied in this work for /I -variation. 

d The strongest H2 component does not coincide with the centre of the HI absorption. 
e Observed with the Hubble Space Telescope. 
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nents; HD is abundantly present and can be included in the analysis. This makes 
J 1237 +064 a target with a high potential. Further details on the system in [ 38 1. 

• Q1443+272 is a special case, and its unusually high redshift (z a bs = 4.22) alters 
some of the arguments. The system has a very low brightness, but that is partly 
compensated because the H2 lines shift further into the visible, where scattering 
is less and CCDs are more sensitive. In addition the high redshift shifts many 
more H2 lines into the observable range. However, for these high redshift ob- 
jects the Lyman-a forest becomes denser, which is a disadvantage. In summary 
this system has a good potential, if sufficient observation time is made available. 
Further details on the system in [39|. 



5 Conclusion 

We have reviewed the status of investigations into variation of the proton-electron 
mass ratio from analyses of highly redshifted H2 absorption systems. As of now 
five results have been obtained putting a tight constraint at the Afi/fj, < 1 x 10~ 5 
level. We present spectroscopic data on the fifth H2 absorption system, Q2348 — 
Oil, and have shown how a complex velocity structure can be unraveled using the 
comprehensive fitting method. An attempt is made to list in full all relevant quasar 
systems where H2 has been observed so far, followed by an assessment as to which 
systems could produce a competitive result on Ajj./ji, provided that high-quality 
well-calibrated observations are performed. 
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stalled at the VLT Kueyen UT2 on Cerro Paranal, Chile. The authors wish to thank John Webb, 
Julian King and Steve Curran (Sydney), Adrian Malec (Melbourne) and Freek van Weerdenburg 
(Amsterdam) for fruitful discussions. 
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